home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / TOEPLZ.DEM < prev    next >
Text File  |  1991-04-29  |  840b  |  41 lines
  1. PROGRAM d2r7(input,output);
  2. (* driver for routine TOEPLZ *)
  3. CONST
  4.    n=5;
  5.    twon=10;    (* twon=2*n *)
  6. TYPE
  7.    glnarray = ARRAY [1..n] OF real;
  8.    gltwon = ARRAY [1..twon] OF real;
  9. VAR
  10.    i,j : integer;
  11.    sum : real;
  12.    r : gltwon;
  13.    x,y : glnarray;
  14.  
  15. (*$I MODFILE.PAS *)
  16. (*$I TOEPLZ.PAS *)
  17.  
  18. BEGIN
  19.    FOR i := 1 to n DO BEGIN
  20.       y[i] := 0.1*i
  21.    END;
  22.    FOR i := 1 to 2*n-1 DO BEGIN
  23.       r[i] := 0.1*i
  24.    END;
  25.    toeplz(r,x,y,n);
  26.    writeln('Solution vector:');
  27.    FOR i := 1 to n DO BEGIN
  28.       writeln('x[':7,i:1,'] :=',x[i]:13)
  29.    END;
  30.    writeln;
  31.    writeln('Test of solution:');
  32.    writeln('mtrx*soln':13,'original':12);
  33.    FOR i := 1 to n DO BEGIN
  34.       sum := 0.0;
  35.       FOR j := 1 to n DO BEGIN
  36.          sum := sum+r[n+i-j]*x[j]
  37.       END;
  38.       writeln(sum:12:4,y[i]:12:4)
  39.    END
  40. END.
  41.